Spin Gap in Chains with Hidden Symmetries 
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We investigate the formation of spin gap in one-dimensional models characterized by the groups 
with hidden dynamical symmetries. A family of two-parametric models of isotropic and anisotropic 
Spin-Rotator Chains characterized by 57/(2) x 57/(2) and SO (2) x SO (2) x Z2 x Z2 symmetries 
is introduced to describe the transition from SU(2) to 50(4) antiferromagnetic Heisenberg chain. 
The excitation spectrum is studied with the use of the Jordan- Wigner transformation generalized 
for 04 algebra and by means of bosonization approach. Hidden discrete symmetries associated with 
invariance under various particle-hole transformations are discussed. We show that the spin gap 
in SRC Hamiltonians is characterized by the scaling dimension 2/3 in contrast to dimension 1 in 
conventional Haldane problem. 
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More than 20 years ago Haldane [l| made a conjec- 
ture that the properties of spin 5 Heisenberg antiferro- 
magnetic (AF) chains are different for integer and half- 
integer spins. Namely, the excitations in the Heisenberg 
AF chains with half-integer spins are gapless, whereas 
for integer spins there is a gap in the spectrum (Hal- 
dane gap). While the first part of Haldane conjecture has 
been proven long time ago (see 0,0); the second part, 
although being confirmed by many numerical Q and ex- 
perimental studies and tested by some approximate 
analytical calculations 

0-E3 

remains a hypothesis. The 
problem of SU(2) Heisenberg chains has been attacked 
by the modern tools as e.g. bosonization la-Ja (see also 
book dlj), various numerical methods [il \l2L H?3 ] and 
recently proposed fermionization by means of Jordan- 
Wigner transformation for higher spins 14] . However, 
the main focus of interests has been put either on 5=1 
chains characterized by SU(2) symmetry or on A^-leg 
ladders described in terms of dynamic SO(N) groups 
(lif . There have also been made several conjectures con- 
cerning spontaneous discrete symmetry breaking in 5=1 
chain models associated wit h e.g . existence of hidden Z2 
and Z2XZ2 symmetries |l3l Il6j. Nevertheless, the gen- 
eral question about the nature of spin gap is still open. 

In this paper we propose yet another approach to the 
spin gap problem. It is based on investigation of a family 
of two-parametric Hamiltonians described by dynamical 
groups 0. This family includes conventional two- leg 
ladder and several models intermediate between the lad- 
der and the chain. Here we concentrate on the most 
instructive example of a " barbed- wire-like" chain with 
spins 1/2 in each site coupled by the ferromagnetic ex- 
change J± within a rung and the AF interaction J» along 
the leg (Fig.l). The model Hamiltonian is 
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FIG. 1: Spin Rotator Chain. 



model to a case where the states on a given rung form a 
triplet/singlet pair. We call the chain shown in Fig. Q]thc 
Spin Rotator Chain (SRC) (in contrast to the spin rotor 
model [1(1 lisj y Unlike earlier attempts to construct the 
representation of 5 = 1 state out of s = 1/2 ingredients 
0>@, we respect in this case the 50(4) symmetry of spin 
manifold on each rung [lflj ]. As a result, the singlet state 
cannot be projected out. Moreover, it plays integral part 
in formation of the spin gap. We show that the hidden 
Z2 symmetries in this model are the intrinsic property of 
the local SO (4) group of spin rotator on the rung, and 
the symmetry breaking due to nonlocal (string) effects 
results in spin gap formation. These special symmetries 
distinguish our model from N > 2-leg ladder models and 
SU(2) chains. In particular we show also that the scaling 
dimension of a spin gap in SRC differs from that in 2-leg 
ladder. 

New variables on a rung are introduced to keep track 
on S = 1 properties . We define Si = Sij + s*2,i ; R% — 
si,i — S2,%, where Si stands for a triplet 5=1 ground 
state and singlet 5 = excited state. The operator R 
describes dynamical triplet/singlet mixing 

[HEI. Then 
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This model is a natural extension of the 5=1 chain 



J- [SiSi+i + SiRi+i + (S^R) 

i 
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where the set of operators Si,Ri fully defines the 04 al- 
gebra in accordance with the commutation relations 

[R?,S?] = i6 ij e af}l Rl (3) 

where e Q ^ 7 is the totally antisymmetric Levi-Civita ten- 
sor and Casimir constraints on each sites are given by 

(5,) 2 + (R t ) 2 = 3, (Si ■ Ri) = 0. (4) 

In order to characterize low-lying excitations in SRC 
we propose a fcrmionization procedure, which extends 
Jordan- Wigner (JW) transformation to 5*0(4) group, 
and a bosonization formalism based on this procedure. 
Our method incorporates JW transformation for 5 = 1 
proposed by Batista and Ortiz (BO) in Q. The BO 
representation is however redundant and requires a con- 
straint overlooked in [l4(. The relationships between 
50(4) JW representation and BO representation is dis- 
cussed below in some detail. 

We begin with a single-rung dimer problem. A two- 
component fermion (a+6+) basis representing 5-operators 
is introduced as follows (5 ± — S x ± iS y ) 

5+ = a++e i7ra+ °6+, 5- = a+be-' lnata , S z = 0+0+6+6-1. 

The complementary representation for R generators is 

R+ = a+ - e Wa 6+, R- = a - 6 e - 47ra+Q , R z = a+a - 6+6. 

This representation satisfies commutation relations © 
for the 50(4) group and preserves Casimir operators 
. The advantage of two- fermion formalism in compar- 
ison with two independent JW transformations for each 
s = 1/2 is that the latter requires an additional Ma- 
jorana fermion to provide commutation of two spins on 
the same rung. Two-component spinless fermions may 
be combined into one spin fermion, which is most conve- 
niently done by the definition 

/ T = (a-&)/V2, fl = (a + b)/V2. (5) 

In order to generalize one-rung representation for a linear 
chain of rungs we introduce a "string" operator Kj 

Kj = exp[j7r 2J »Vfc] = JJ(1 - 2n Tfe )(l - 2n lk ), (6) 

k<j,a k<j 

(n>cr — fifa)- As a result of JW transformation the 
50(4) generators acquire the following form 

Sf = V2 (4(1 - nii )Ki + K}f u (l - ny)) , 

Sr = (5+) f , S] = n rj -n u , (7) 

R t = ^{fh n uKi + K]f u n h ), 

l\, Cry- R]=fyh+flifrr (8) 



Part of the representation J7J) describing 5=1 coin- 
cides with BO representation. Nevertheless, since 5 2 
is no more a conserved quantity, being defined by 
5?=2[1— nyriij], the projection of 50(4) group on 5=1 
representation of SU(2) group requires an additional 
Hubbard-like interaction responsible for the hidden con- 
straint overlooked in BO paper JJ^j. When the 5=1 sec- 
tor is fixed, three states (n-f,nj.), namely (1,0), (0,0) 
and (0,1) determine three-fold degenerate triplet state 
whereas the doubly occupied state (1,1) stands for a sin- 
glet separated from the ground state by the gap A=J±. 
The Hamiltonian (01 is fcrmionized by means of purely 
ID string operator Kj © in contrast to meande ring 
strings proposed for the theory of 2-leg ladders (see |2(| 
and references therein). 

The Hamiltonian of anisotropic XXZ SRC model is 
H = H \\ +EiH±,i, where 

H \\ = -j- E [#^1 + S ? R ^ + (s~r)+ h. c ] 

i 
i 

H± ti = § (R+K + Ri R t) + J -f{Rt) 2 - (Ri - 

There exists a set of discrete transformations keeping the 
Hamiltonians |J2J and @ intact and preserving commuta- 
tion relations © and Casimir operators QJ. In general, 
these transformations are described by the matrix of fi- 
nite rotations characterized by Euler angles 9, tp, <p, tp for 
the case of SU{2) x SU{2) or 50(2) x 50(2) x Z 2 x Z 2 
groups. An example of such transformation is 

5+ -> R+, S~ -» R-, S z -» S z , R z -» R z . (10) 

being a t/(l) x £7(1) rotation in " S x ~ and "5 y -i? y " 
subspaces. This is in fact a particle-hole flavor transfor- 
mation fi — ► /|, /| — > /J. On the other hand, it corre- 
sponds to replacement b — ► — & thus manifesting hidden 
Z2 symmetry. This means that an additional gauge fac- 
tor exp(i#) with 9 = ±7r appears in a fermion operator 
characterizing " free ends" of rungs in SRC chain. Other 
examples are (/ T -> / x ) and (/ T -> /+,/; -> /+). The 
latter one corresponds to a particle-hole transformation 
(a — ► a+, 6 — ► 6+) in the non- rotated fermion basis. 

After JW transformation in a — b basis ©-© the 
Hamiltonian is written as follows 

H\\ = J\\ ^2 ( a ! fl i+i + a l+i a ») cos(7rn-) 

i 

+ j iiE(<-£) (^1-5) ( n ) 

and H± — J^. j with 
#M = + bUi) ~ Jl{< - J)(nf - h (12) 
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FIG. 2: Dispersion law for hybridized spin fermions c±. 

where the shorthand notations n a = a^a, n b = b% and 
cos(7rn b ) = Reexp(±i7rn b ) = 1 — 2n b are used. Below we 
consider the domain J± <C J\\ where strongest deviations 
from conventional Haldane gap regime 0,0 are antici- 
pated. In the limit J± = our SRC model reduces to 
an s = 1/2 AF chain, the gauge factor cos (7m ) = ±1 
is a fictitious random variable which can be eliminated 
by S x — ► — S x and S y — > — S y on the corresponding site. 
This situation is similar to the so called Mattis disorder 
|2lj | where the randomness in interaction is removed by 
proper redefinition of spin variables. 

The kinematic factor ~ cos(7rn|) in H x Ijllj) can be 

eliminated by a unitary transformation H=U^ HU with 
U=exp(iTr^2 l j>i n j n i )• Then H\ and Hi remain un- 
changed and term acquires the string form 

Oft = [alhe-^^M^M + h.c)j . (13) 

The s = 1/2 chain is represented in terms of a half-filled 
band of fermions. Since interactions I jlljl and (j!2jl do not 
change the occupation numbers for each color, we expect 
that interacting case is also represented by two half-filled 
bands (see below). We note that the Hamiltonian H\\ in 
@ possesses U(\) x £7(1) symmetry whereas only one 
local U(l) associated with b- fermions exists in (llljl due 
to non-local character of JW transformation. 

Let us consider the XY\\ — XY±_ model (J,| = = 0). 
We spilt the first term in 1111) into the bare hoping and 
the kinematic term ~ Jifn^af cii+i + h.c) playing part of 

effective interaction H^J . One gets after diagonalization 
of the hopping term 

Ha= £ a(p)4 iP c AiP (14) 

p,A=± 

with c + = u + a + c_ = u + b — w_a, 

4(p) = ±s ± (p)/(e + (p)-e.{p)), (15) 

e±(p) = J\ ] cos P ±[{J x cosp) 2 + {Jl) 2 ] 1 ' 2 . (16) 

The chemical potential ~q=0 is pinned in the gap. Thus, 
the mixing term fixes global phase difference for a — b - 
fields. The remaining symmetry is local Zi^Z^- 



We represent H^J in terms of new variables c± by 
expanding the Hamiltonian (jlljl in the vicinity of two 
Fermi points of non-hybridized system 

H?J = \ ^(?)P^'.«(?)A^,a'(-g) (17) 

where the operator p^i is given by 

P^',a(q) =Y C a,V,k-q/2 C <*,ti',k+q/2- (18) 
k 

Its diagonal part is the quasiparticle density. The oper- 
ator k vv i is defined as 

Ku',<x{q) =~a22 k ^ a .v,k-q/2 C ^',k+q/2, (19) 
k 

while its diagonal part is K vv — divj,,„ = —dtp vv . 

In expressions Ijl8j) . Ijl9jl the index a = ± stands for 
"old" Fermi surface points kp = ±7r/2 (we take a unit 
lattice spacing), k is measured from kp. Indices fx, p', 
v, i/'=± denote the branch of fermions c±. We used the 
property of u± iCt (±7r/2)wl/v^. The tensor ^ for these 
scattering processes has the form 

C' = ^(V + ^OOW-^V)- (20) 

We analyze Ijl7jl in terms of <?-ology approach j^l 
classifying various terms in g^i(q) as forward, back- 
ward scattering and Umklapp processes. We see, first, 
that if |g|<C7r/2, and g~±J x } both diagonal and off- 
diagonal matrix elements of h. vv i vanish in accordance 
with Adler's principle j2^|. Thus, the forward scattering 
processes leading to small renormalization of the coupling 
r-j (J±) 2 /J x are irrelevant. The backward scattering pro- 
cesses (±7r/2^=F7r/2) result in a reduction J^^J^ 1 J 1 ^ 
of the effective coupling (0<7<1 is a constant). To get 
this estimate we cut logarithmic corrections to the cou- 
pling constant by A m i„~( J]_) 2 / J\\ where A m .;„ deter- 
mines the gap in the density of spin-fermion states e±. 
However there is yet another energy scale A~ associ- 
ated with the gap in a two-point particle-hole correlator 
with zero total momentum of the pair. This energy scale 
is attributed to the gap separating 5=0 excited state on 
a rung from the triplet state. The crossover between two 
energy scales will be discussed elsewhere. The Hamil- 
tonian Ijl7jl allows also "inter-band" Umklapp processes 
determined by the off-diagonal elements of p^i and K vv > 
and responsible for periodicity Q=2ir. These processes, 
associated with the transfer of pair of quasiparticles over 
the gap do not change the leading term in Ijl6jl . 

The above arguments are complemented by the 
bosonization calculations for the strongly asymmetric 2- 
leg ladder with finite Fermi velocity Ub in the b subsystem 
which may be turned to zero in the end of scaling pro- 
cedure. The continuum representation for spin operators 
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si,s 2 in (JTJ reads [UHl] (we denote (t=o(l), 6(2))) 

sf(x) ~ e ±4e! (cos(7ra;) +003(20,)), 

s*(x) — TT^ 1 d x <p l + cos(ttx + 2(j>i) (21) 

with canonically conjugated variables fa and IL„=9 x 0j. 
Keeping only most relevant terms in the rung interaction, 
J", we arrive at the conventional equations of Abelian 
bosonization for the spin Hamiltonian Q 

+ J* cos(# a - 6 ) + J* cos 20 a cos 20b] (22) 

with if =1/2 and J±^.Ub^u a =irJ\\ /2 for J^J^^I. 

To find the scaling dimension of the gap we start with 
the case J^=0, J±=J±^=0. Using the scaling proce- 
dure [x — > Ax, t — > At), one has jj_ — > J±A 2 ~@ where 
/3/2=K/2 is a scaling dimension of cos20i. The renor- 
malization of 6-component stops when the renormalized 
J± becomes comparable with the lower scale of the en- 
ergy Ub- The corresponding scale A = defines the 
first correlation length £5 = 'J±) l '^ 2 ~^ and the first 
energy gap = Ub^ 1 — Ub(J±/Ub) 1 ^ 2 ~ 13 ^ ■ At the sec- 
ond stage of the renormalization, with frozen (cos 2<pb) ~ 
£ b the factor cos 2(f> a undergoes further enhance- 
ment. The procedure halts when the renormalized 
amplitude Jj_ compares with u a at JiA 2_,3 / 2 £ b ~ 
u a , which defines a second correlation length £ a = 
(^ /2 Ua/J±) 1/(2_/3/2) and the second gap A a = u a £-\ 
In our particular case (3=1 these formulas simplify as 
follows: 

£,b = (u b /J±), A fc = J± 

£ Q = &(«o/«6) 2/3 , A a = J±(u a /u b ) 1/3 . (23) 

One may decrease Ub in the regime of frozen down to 
Ub^Jj_- Then A a ~ J11 ( Jj_/ J11 ) 2 / 3 . Further decrease of 
does not change the exponent 2/3 of the spin gap fully 
determined by the scattering on the random potential 
cos 2<fi a [2^ . The two-stage renormalization procedure is 
essential for understanding the SRC model. In the limit 
Ua ~ Ub, Eq. lll'-'il) leads to standard scaling of the spin 
gap A ~ J± (see e.g. H^). 

In the case Jf^O, J\=0 the scaling behavior of the 
spin gap A~ J11 {J±/J\\ ) 2 ^ 3 is determined by the backward 
scattering processes of the field a on the random potential 
associated with fluctuations of cos# Q . 

The fully isotropic case, J±=J±=Jx, might be ex- 
pected to yield the same estimate A~J,V 3 (,7j_) 2 / 3 . A 
refined analysis (see, e.g. [2?|]) including the less relevant 
terms in (I17|) may correct the gap values, but does not 
change this estimate. 

To summarize, we introduced a new ID model inter- 
mediate between the spin 5=1 chain and the 2- leg lad- 
der. Our SRC possesses special hidden Z2 symmetries 



connected with discrete transformations in a 6D space of 
50(4) group characterizing the spin rotator. The SRC 
chain is mapped on the two-component unconstrained 
interacting fermions by means of 04 JW transformation. 
Two fermion fields are characterized by sharply different 
dynamics. One of two fields is frozen at /c^±7r/2 and the 
scaling dimension (3 of the rung operator exchange J± is 
(3=1/2 instead of (3=1 [2^, (2^| in a conventional Haldane 
problem. As a result, new scaling "two third" law for the 
spin gap arises. 
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